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The screw instability of magnetic field is a mechanism for prohibiting a generation of strongly 
twisted field lines in large scales. If it can work in black hole magnetospheres, the global axisymmetric 
structure and the main process of energy release will be significantly infiuenced. In this paper, 
we study the instability condition in the framework of the variational principle, paying special 
attention to a stabilizing effect due to field-line rotation against the screw-unstable modes satisfying 
the well-known Kruskal-Shafranov criterion. The basic formulation for the stability analysis of 
rotating force-free fields in Kerr geometry is provided. Then, for the practical use of the analytical 
method, the stationary configuration is assumed to be cylindrical, and we treat the long-wave mode 
perturbations excited at large distances from the central black hole, where the strength of gravity 
becomes negligibly small. This allows us to derive clearly the new criterion dependent on the 
rotational angular velocity of magnetic field lines, and the implications of the results obtained for 
magnetic activities of jets driven by a black hole are briefly discussed. 

PACS number(s): 04.70.-s, 95.30.Qd, 98.62. Js 



I. INTRODUCTION 

The global structure of black hole magnetospheres involving axisymmetric magnetic field and plasma injected from 
an accretion disk has been extensively investigated for explaining various observational features of active galactic 
nuclei (AGNs) (see, e.g., 0] for a review). The magnetic field lines which are assumed to be frozen to the rotating 
plasma may be twisted and rotating, which means that the components of toroidal magnetic field (i.e., poloidal electric 
currents) and poloidal electric field are generated. The magnetic tension of the toroidal field can play an important 
role of forming cylindrically coUimated jets over significantly long distances, and if the magnetic field lines thread 
a rotating black hole in the central region of AGN, the so-called Blandford-Znajek mechanism |^,^ is believed to 
become a possible process for powering the plasma jets through Poynting flux. 

However, it is also well-known from the field of plasma physics that the magnetic field configurations with both 
poloidal component and toroidal one can be screw- unstable [Q|5j, and this process of a current driven instability 
makes the magnetic field lines untwist on short time-scales and causes a sudden release of magnetic free-energy to 
plasma's kinetic energy. In conventional plasma confinements, the screw instability is claimed to work for long-wave 
mode perturbations, if the equilibrium magnetic field has a toroidal component exceeding the limit given by the 
Kruskal-Shafranov criterion. Then, according to this instability condition one may discuss that the power of the 
Blandford-Znajek mechanism which is crucially dependent on the strength of of toroidal magnetic field is significantly 
suppressed and the mechanism for acceleration and coUimation of jets also becomes less efficient. Further, a 
flare-like energy radiation is expected to occur in the screw- unstable magnetospheres . After a part of the magnetic 
energy is released through the process of the screw instability, the magnetic field lines threading a rotating black hole 
may be twisted again for generating the toroidal component stronger than the Kruskal-Shafranov threshold, and the 
fiaring may be quasi-periodically observed. 

Thus, the screw instability is an astrophysically interesting process which can significantly infiuence the global 
axisymmetric structure and the main process of energy release in black hole magnetospheres. For the application to 
the astrophysical problems, however, the instability condition should be investigated in more detail. In particular, 
an essential feature of black hole magnetospheres is that the twisted magnetic field lines are rotating. In this paper, 
our main purpose is to consider the effect of the field-line rotation in the stability analysis of stationary magnetic 
configurations and to derive the new criterion of the screw instability. In Sec. II, in the framework of the force-free 
approximation in Kerr geometry, we give the basic equations for the axisymmetric stationary fields and the non- 
axisymmetric linear perturbations, following the usual procedure to derive the variational principle For the 

practical use of the analytical method, in Sec. Ill, the stationary configuration is assumed to be cylindrical, and we 
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treat the long-wave mode perturbations excited at large distances from the central black hole, where the strength of 
gravity becomes negligibly small. Then we can clearly see that the field-line rotation has a stabilizing effect against 
the screw-unstable modes satisfying the Kruskal-Shafranov criterion. This allows us to reconsider the conclusions 
based on the Kruskal-Shafranov criterion as a first step toward more complete investigations using the fully general 
relativistic formulae given in Sec. II. While any direct influences of Kerr geometry are missed in the new criterion 
obtained here (which is relevant to force- free rotating magnetospheres far distant from various central objects), it is 
possible to discuss the implications for magnetic activities in black hole magnetospheres, by estimating the toroidal 
component of the stationary magnetic field under the assumption that it is generated as a result of the black hole 
rotation. Section IV is devoted to a brief comment on this problem. Throughout this paper, we use the geometric 
units with G = c ~ 1. 



II. STABILITY ANALYSIS IN KERR GEOMETRY 

In the framework of the 3-1-1 formalism in Kerr geometry we study a current driven instability of electromagnetic 
fields under the force-free approximation which will be valid in black hole magnetospheres except near the event 
horizon and the (inner and outer) light cylinders |^. The Kerr metric is written in the Boyer-Lindquist coordinates 
as follows, 



where 





-a^dt^ + h,j{dx^ - P^dt){dx' - P^dt), 


(2.1) 
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p' 


^r^ cos^ 9, A = + - 2Afr, 


(2.4) 


Y? = 


{r'^ +a'^ f -a'^Asin^e, to sin 6*, 
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(2.5) 



and M and a denote the mass and angular momentum per unit mass of the black hole, respectively. Then, the 
Maxwell equations and the force-free condition for zero-angular-momentum observers (ZAMOs) are given by 

{dt+Lfi)B^-VxaE, (2.6) 
(dt + Lp)E ^VxaB- aj, (2.7) 
pE + jxB^O, (2.8) 

where denotes the Lie derivative along /3*, B and E are the magnetic and electric fields, p and j are the electric 
charge and the current densities multiplied by 47r ||^. 

If the current and charge densities are eliminated from the Maxwell equations, we have the equation 

dt{B X E) + B X LpE + (LpB) x E + a(V • E)E = B x {V x aB) + E x {V x aE). (2.9) 

Because the force- free fields may be regarded as a magnetically dominated limit in magnetohydrodynamics (MHD), 
we introduce the plasma velocity defined by 

E + vxB^O, (2.10) 

and we have the frozen-in law written by 

{dt+Lp)B = ^ x{vxaB). (2.11) 

These are the basic set of dynamical equations for the force- free fields. 

For axisymmetric stationary configurations of electromagnetic and velocity fields denoted by Bo, Eq and wq, we 
obtain the well-known relations. For example, the poloidal component B^ of the magnetic field is written by the 
magnetic stream function 5* as follows. 
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27r-!i7 ' 



(2.12) 



where eV = [1 / 'uj){d / d(j)) means the toroidal basis vector which has unit norm. In stationary rotating magnetospheres 
the velocity field is given by 



(2.13) 



where Vlp = ^f{'^) is constant along a magnetic field line of '^{r,6) = constant, and we call it the angular velocity 
of a magnetic field line. This rotational motion generates the poloidal electric field Eq = —vq x Bq. Further, if the 
magnetic field lines are twisted, the poloidal current / — /(^') (which is chosen to be positive for downward flows) is 
generated according to the equation 



(2.14) 



where B^ is the toroidal magnetic component. Then, Eq. (|2.9D reduces to the Grad-Shafranov equation for ^ ||l 



Utt^ ^ dl 



1^=0, 



(2.15) 



which should be solved by giving the functions /(vl/) andJ7p(^) from boundary conditions at a plasma source. Though 
we do not try to solve the highly nonlinear equation (2.15) for background fields, it will be used for simplifying the 
equations for linear perturbations in the stability analysis in the following. 

Now, we consider small pert ur batio ns Bi, Ei and vi from the stationary configurations. It is easy to derive the 
linearized versions of Eqs. (U), ( px| ) and ( pi] ) as follows. 



El = -Wo X Bi - wi X Bo, 
{dt + Li3)Bi = V X (wi X aBf) + vq x aBi), 
Bo X dtEi +BoX LpEi + B^ x L^Eq + a [(V • Eo)Ei + (V • Ei)Eq 

= i?o X (V X aBi) + i?i X (V X ai?o). 



(2.16) 
(2.17) 

(2.18) 



Because the background fields and geometry are stationary and axisymmetric, the perturbed fields can have the form 
of exp(iTO(/) — icri), where m is an arbitrary integer and a is an arbitrary complex number. The key point for analyzing 
the linear equations is to introduce the plasma displacement ^ defined by vi — ^icr^, for which we can give ^ = ^ x Bo- 
Then, from Eqs. ( 2.16| ) and (2.17), the perturbed fields Bi and Ei are explicitly solved as follows, 

Bi 



v{\7 xaC)+ ;^weT, 



E^ =v 



■w{U,F — i^) 



a — m^F 



(V X aQ) X ct 



(2.19) 
(2.20) 



where v — a jia — mflp). 

For mathematical simplicity in later calculations, we assume that is constant all over the field lines, i.e., the 
second term of the right-hand side of Eq. ( 2.19| ) vanishes. Then, by substituting Bi and Ei into Eq. (2.18), our 
problem reduces to the eigenvalue problem for the cigcnfrequency cr and the eigenvector C (or ^ ) , which has the form 



n=0 



(2.21) 



Because we use the variational principle to consider the eigenvalue problem, here we express the differential operators 
An in terms of the integrals of arbitrary displacement vectors ^ and fj as follows. 



(2.22) 



where stars and h denote complex conjugates and the determinant of the metric hij, respectively. Some parts of 
the integrands ff * ■ (A„^ ) reduce to the total divergence written by V • 5, which is assumed to vanish owing to the 
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condition such that no perturbations are generated at the boundaries. Then, through the lengthy calculations using 
C and X = 77 X i?o instead of and 77, we arrive at the results 



B, 







-{rtF-Lo) X * • V(n7aCT) - C • V(n7axT) 
a L 

^2 . 



(2.23) 



(2.24) 



bo — m^w^x * ■ C + 



-V(n7Q!XT) ■ V(ci7aCT) 



I '-Bo • (x X C ) (f7i? - u) 



V{waCT) - C • V(roaxT) 



(V X ax *) • (V X aC ) - 2^ ( — ) X * • (V X aC ) + C • (V X ax *) 



(2.25) 



where C,t and xt are the toroidal components of Q and x measured in the orthonormal bases. These equations clearly 
show that the operators An are Hermitian. 

If the vector x is chosen to be equal to the coefficients a„ appearing in the dispersion relation 



a2(J + aiCF + ao = 



(2.26) 



become real numbers. Because a2 is positive definite, it is clear that the sufficient and necessary condition for the 
screw instability of stationary configurations is 



W = Wp + Wb,<Q, 



(2.27) 



. Even if Wp is not positive 
yy virtue of the existence of the 



where Wp = —ao is the so-called potential energy for the small displacement ^ Q 
semi-definite, stationary configurations in rotating magnetospheres can be stabilized 
rotational term Wp. = a\/Aa2. 

It should be remarked that a„ are the functional of C which is perpendicular to Bq. This motivates us to introduce 
the poloidal orthonormal basis vectors defined by 
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2tiw\BI\ ' 



(2.28) 



in addition to the toroidal basis vector eV with the condition e|| — y. ct- Then, decomposing the vector C into 

C = C||e|| +aei + CTeT, (2.29) 

we obtain 

Ct = -UBl Cll = e±So^ (2.30) 

for Bq = -BgCii -|- BqEt and £,± = e± ■ ^. Hence, we can treat the two functions and as the trial functions to 
calculate the eigenvalue a. 

We have established the basic formulae to use the varia tional principle for the stability analysis in Kerr geometry. 
The next step will be to apply the method to Eq. ( 2.26| ) with the purpose of determining the trial functions and 
obtaining the instability criterion expressed by stationary background fields. Unfortunately, this becomes a very 
complicated task owing to the geometrical terms involved in a„. Then, our interest in the next section is focused on a 
stabilizing effect of flp in rotating magnetospheres, considering simple stationary fields and the perturbations given 
in flat spacetime. 



III. STABILIZING EFFECT IN ROTATING MAGNETOSPHERES 



As a mechanism of jet formation the existence of collimated field lines threading a black hole is astrophysically 
interesting. The jets would propagate along the poloidal magnetic field lines to large distances, keeping the collimated 
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structure under an action of the toroidal magnetic field. Here, the stationary magnetic field Bg is assumed to be 
cylindrical as a typical configuration of coUimated field lines. We discuss the instability of such a global structure due 
to the perturbations ^ with a wave length (in the vertical direction along the filed line) much longer than the horizon 
scale. (Long- wave modes are expected to become dominant in the screw instability H].) Then, the integrals in the 
vertical direction appearing in a„ should extend to a very large scale such that L ^ M , where gravity becomes very 
weak. This allows us to estimate approximately the integrals a„ by giving fields in flat spacetime. 

Because of the cylindrical configuration of the stationary magnetic field, it is convenient to use the cylindrical 
coordinates denoted by vc, (j) and z, for which the Euclidean metric hij is given by 

hijdx^dx-' = dvo"^ + w^d^i^ + dz^ , (3.1) 

and the orthonormal basis vectors previously introduced become 

6*11 = 6*2, e_i = e^. (3.2) 

The magnetic components are assumed to be dependent on -w only as follows, 

Bo = B{w)e^ + U{w)eT, (3.3) 



and for a = 1 and co — the Grad-Shafranov equation ( 2.15 ) reduces to 



(1 - zD^nj,)BB' + UU' = 2Tun%B'^ - —U^, (3.4) 



where primes denote the differentiation with respect to w. To derive Eq. (3.4), we have also used the relation 

dl _ {vjUy 
d^ AnzuB 



(3.5) 



for the poloidal current /. If the toroidal component U is given, we can determine the distribution of the poloidal 
component B. The simplest example giving a constant B would heU = izuflpB. 
Now the perturbations with the eigenfrequency a can have the form 

— — exp(im(/< + ikz — iat), 

C± ~ Cxa^ — C('^) exp(«m0 + ikz — icrt), (3.6) 

where k is the wave number in the z-direction, and the minimal value of k becomes approximately L^^ in the 
magnetosphere with the scale L. To obtain a relation between the two trial functions ^, ^, we consider the variation 



of Eq. (2.26). As will be seen later, the variational principle becomes consistent if fj is chosen to be an eigenvector 



with the eigenvalue a* (which may be complex), i.e. 



Vi- = = ^(^) exp(ir7i(/) + ikz — ia*t), 

= Xro = x(^) exp{im(j) + ikz - ia*t), (3.7) 



where 77 — ^(nj, a*) and x — C('^7 



Using these eigenvectors ( and x, the coefficients a„ in Eq. (2.26) reduce to 



a„ = 2ttL J wdw bn- (3.8) 
with L as the length in the z-direction. It is easy to calculate the integrands 62 and bi with the results 

b2^X*C+iB' + U^)v*i, (3.9) 

61 = iflp [Ci^B?]'* - Br]*) - x*{^B^' - B^) + 2i{mB - kwU)BifS,] . (3.10) 

To calculate 601 the partial integration of the terms containing (?7*0' becomes necessary. By omitting the contribution 
from the boundary, we have 

b. = Pr/r + Qrn + !!^±^ Ux* - U + , (3.11) 
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where 



P 

Q 

Si- 



{mU + kvjBf 
{mU + kwBf 



1 



0^% 



(1 - - k^zu'^)B^ - 



-1 - 



2k^k'vj'B' -m'U^) 



2 ,2 2^o2 2k^w^imBB' + 2B^) 



{mB - kwU)vj£_' + {mB + kwU)^, 
(mB — k'cuU)wri' + (mB + kzuU)?]. 



(3.12) 



(3.13) 

(3.14) 
(3.15) 



I t sho uld be noted that no derivatives of ( and x appears in the integrands. Then, taking the variation 6 of Eq. 
( 2.26 ) and putting = S-q* = 0, we require that the coefScients of 5C, and 5x* vanish, because they are algebraicaUy 
independent and arbitrary. This procedure leads to the results 



and 



(m^ + fc^tn^ - cr^ro^)^* = -iS"* + inazu^i-cuBri'* - Sr;*), 



while the latter is automatically satisfied by virtue of 77 = ^(ci7, a*) and x — Ci^_,<^*)- 

Since ^ has been fixed by the variational method, the dispersion relation (2.26) is now given by 



/ 



wdw{FTj'*C + Gtj*^ + H„) = 0, 



where 



F = - 



{mU + kvjBf 



(3.16) 
(3.17) 

(3.18) 
(3.19) 



G 



{mU + kvjBf 



and the term dependent on a is 



H„ = 



m 



k\ 



m'^niiw^B^y 



(3.20) 



l2 + k^Vj"^ 



k^w^ 



npBiwri'* - 7]* 



k'^w'' 



- flpBizuC - 



[(cr - mVLpfB'^ + {<tU + kru^FB)^] ri*£_. 



(3.21) 



If a is real, the term Ft]'*^' is negative semi-definite. Then, if it becomes dominant in Eq. (3.18), the third term 
can be positive for giving consistently a real eigenvalue. This means that the most dangerous mode to make a 
complex is obtained by choosing ^ in such a way that the integral involving the term Frj'*S/ vanishes. In fact, for 
f2i? = R], the screw instability has been found to occur for the mode given by 



1 , for -uj < wq, 

{roo(l + fi) - w}/ (wqh) , for njQ < in < nJo{l + fi), 
, for ti7o(l + /i) < ro. 



(3.22) 



where wq is the first zero of mil + kwB which is assumed to be present in the range < n7 < 00. Taking the limit 
/i ^ 0, we can easily estimate the integral of Fri'*£,' to be of the order of ^. (This trial function is define d without 
including cr, and so we have ^ = 77.) We use the same trial function to analyze the dispersion relation (3.18) extended 
to the case of flp ^ 0. 

Let us assume that B is everywhere positive. For regularity of the toroidal field on the symmetry axis, U should 
approach zero at least in proportion to m in the limit zu ~* 0. Further let the wave number k be positive. Then, the 
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sign of m is chosen according to the requirement such that mU < 0. By virtue of the steep gradient of £_at nj = -toq, 
the integral of containing the derivatives becomes of the order of in contradiction to Eq. ( 3.18 ), unless 
the eigenvalue a is suitably determined. It is easy to see that such a divergent term vanishes, if a is given by 



a — mf2 + (Ti . 

The small correction cti turns out to be of the order of /i^^^, because in the limit /i 



(3.23) 



Eq. (3.18) reduces to 



-J, 



(3.24) 



where ~ B{zuq), and we have 



J = {mU + kmBf 



1 - r^i 



1 



m^(l — 



k^VJ^ 



,212 ■ 



(3.25) 



Though the eigenvalue of ai is no t determined unless the gradient /i^^ of ^ at ti7 = vjq is precisely given, it is easy to 
see the sign of erf from Eq. (3.24). We note that the stability of stationary fields may drastically change if the light 
cylinder (i.e., ^Lpw = 1) is present in the range < ct7 < vjq. However, our simple force-free model would not remain 
approximately valid in a region beyond the light cylinder where plasma inertia becomes very important. Hence, we 
would like to limit our stability analysis to the fields inside the light cylinder, i.e., vjq^f < 1, and we can conclude 
that the stationary fields given by B and U become unstable if the condition 



dw 



-J < 



(3.26) 



is satisfied. This is equivalent to Eq^ ( 2.27 ), if W is estimated by using the eigenfrequency a = Vlp and the trial 

functions given by Eqs. (3.16) and (3.22). 

From the first term proportional to [mU -{-kzuB)'^ in J the screw mode with |m| = 1 (and mU = — |J7|) turns out to 
be most important for the instability. To compare clearly the instability condition with the usual Kruskal-Shafranov 
criterion, let both k and J^f be positive and small, i.e., kwQ <ti 1 and flpTj^o ^ 1- For the long- wave screw mode 
giving mU — —\U\ in slowly rotating magnetospheres we have 



J -{\U\ - kzuB){k^ + nj,)zu^ 



\U\-kzuB- 



A{k^ 



fc2. 



^'kvuB 



(3.27) 



Then, if ilp < k, it is clear that the requirement \U\ > kwB in the range Q < w < wq makes J negative. Because 
the minimal value of the wave number k is of the order of L~^, we arrive at the Kruskal-Shafranov criterion 



ZD L 



(3.28) 



as the instability condition for a twisted magnetic tube with the width w. The important point obtained here is that 
the minimal value of k is also restricted by Qp, and we can only claim that the screw instability surely occurs if the 
condition 



\U\ 



> ripB 



(3.29) 



is also satisfied. 

If the condition (3.2£) breaks down, the stability becomes a more subtle problem depending on the distribution of 
B and U as functions of vj. For example, let us consider a stationary model such that the poloidal field is nearly 
uniform {B ~ Bq), while the weak toroidal field (i.e., \U\ <C B) is given by \U\ = HqvuBq in the range < w < Wi 
and \U\ — rioti7i_Bo in the range wi < w. Though the weak dependence of B on tu may be calculated from Eq. (^), 
it is not important when we estimate the integral of J. The first zero of |?7| — kvjB appears &t w = wq = floWi/k, 
and the condition \U\ > kwB in the range < vj < wq leads to k < \flo\. Then, through the calculation of the 
integral of J, we can easily show that the stationary model with \ flo\ < flp becomes stable even for the screw mode 
satisfying the Kruskal-Shafranov criterion. 
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However, one may consider a different model with the nearly uniform poloidal field and the weak toroidal field given 

1 /2 

by \U\ — ^IqwBq [l — {w/w2)'^] in the range < n7 < n72 and \U\ = in the range 1372 < vd. Then, if the first zero 
wq is chosen to be very close to 1372, we have |?7(tz7)| ^ |[7(n7o)| = kw^B > kwB (i.e., J ~ —U^Vl\w^ < 0) in the 
range Q < vj < wq, except near zu — tuq. The stabilizing effect oi flp becomes very weak for the stationary model 
with \U\ which rapidly decreases as w increases. 

The conclusion that the latter model with a decreasing \U\ is screw- unstable will be premature when the toroidal 
field vanishes at = ci72 very close to the light cyhnder (i.e., flpW2 — !)• This is because we must use Eq. ( 3.25| ) for 



J, when k is chosen to be very small according to the requirement wq ~ n72. Now we obtain the approximate form 

except near w = vjq, and it is easy to see that the integral of J becomes positive by virtue of Q.f'coq ~ 1. Thus, we 
can expect that the rotating field is well stabilized against the screw mode with k smaller than Vlp, unless a significant 
decrease of the toroidal field occurs in an inner region distant from the light cylinder. 



In summary, we have presented the new criterion (3.26), in whic h J is given by Eq.(p.27|) for long- wave screw modes 



in slowly rotating magnetospheres. If the more precise form ( 3.25 ) of J is used, we can show a strong stabilizing effect 
due to the existence of the light cylinder. It is sure that the Kruskal-Shafranov criterion ( 3.2^ ) is not a sufficient 
condition for the screw instability if the twisted magnetic tube is rotating. Importantly, in rotating magnetospheres 
with Q,p > L^^, the toroidal magnetic field may be stably amplified to the maximum value given by \U\ — flpmB, 



instead of \U\ — wB/L. So, we can say that the stationary magnetic field (3.2) is quite a reasonable assumption 



as a global field structure for the jets except loca l structures as a knot. If the deviation from the cylindrical field is 



serious, the criterion based on the integral (3.26) would not be simply applicable. However, by rough estimations, 
we can expect that the stability effect hy VLp remains valid in modified criterion. This will be crucial for discussing 
highly energetic phenomena in rotating magnetospheres. 

IV. DISCUSSION 

Now we would like to give a brief comment on the application of the results obtained in the previous section to 
black hole magnetospheres. Let us assume that radial magnetic field lines thread a Kerr black hole which is slowly 
rotating with the angular velocity uih- The split monopole magnetic field is given by 

vl^ = -27rco cos6l, B^ ^ cq{VLp - uh) 9 (4.1) 

with a constant co, from which we have 

Bl ~ {np - ujh)wBI (4.2) 

in the region distant from the event horizon. As r increases, the radial magnetic field lines would tend to be collimated 
toward a cylindrical configuration, conserving the magnetic flux w^Bq and the electric current zuBq . Then, the 
stationary cylindrical field, for which the instability condition was discussed, would satisfy the relation 

U c^{np-ujH)vDB, (4.3) 

for the field lines threading the central black hole, i.e., in a limited range such that < m < -cui. 

The constraint on the angular velocity flp oi the magnetic field lines has been studied in |^ using the Kruskal- 
Shafranov criterion. However, if the length L of t he magnetic tube is larger than ^Ip^ , the sufficient condition of the 



screw instability should be replaced by Eq. (3.29), which leads to 

\np~LUH\>nF, (4.4) 

and we obtain 

np<uJH/2. (4.5) 

If the magnetic field rotates too slowly (see, e.g., for such a black- hole driven wind model), the toroidal field 
generated by the black hole becomes too strong to be stably sustained against the screw modes. This may yield 
a flaring event in plasma jets, because the stored magnetic energy is released to the kinetic energy of surrounding 
plasma [0. 
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If the magnetic field lines threading the black hole satisfy the condition 



fiF>WH/2, (4.6) 

a stable configuration becomes possible. In particular, the Blandford-Znajek mechanism |^ can work as a steady 
process of energy extraction from the rotating black hole, ii ftp is in the range ujh/'^ < ftp < ^h- Interestingly, the 
magnetic field can be marginally stable, when the power becomes maximum, i.e., flp = ujh/2. Thus, the jets can 
be significantly powered by the Blandford-Znajek mechanism at large distances L far from the central black hole, 
contrary to the claim in jql. 

The force-free model ([d[) is based on the boundary condition at the event horizon, which may not be justified if 



the effect of plasma accretion is taken into account. If we relax the constraint at the event horizon, the relation (4.3) 
would not always hold. Though the condition Qp — LOpijl of the marginal stability could be modified, it is plausible 
that the ratio ojh /^f is a crucial parameter for determining a dominant energy source which powers jets produced in 
black hole magnetospheres: The screw instability will yield a helical structure of plasma jets and flare-like emission 
through the processes such as shock waves, while the Blandford-Znajek mechanism will provide steady Poynting flux 
for plasma acceleration. The physical mechanism to determine uniquely the parameter loh/^p is also controversial 
If various values of the ratio uih/^^f are permissible for rotating black hole magnetospheres in AGNs, such two 
different types of jets may be observable. 

There exists also the Newtonian solution of Eq. ( ^.15 ) for a rotating, paraboloidal and disk-connected field , 
satisfying the condition Bq — —^IpzuBq at large distances. If the same instability condition holds even for the 
paraboloidal field, such a stationary model would be marginally stable. In this paper, we have considered magnetic 
field lines connecting a central black hole with remote astrophysical loads, in relation to jet formation. As was 
suggested in Q , closed field lines connecting the black hole with a surrounding disk may also become screw- unstable 
(see for a non-rotating model). Thus, the analysis of the screw instability of various stationary fields including 
the effect of strong gravity is a necessary task for understanding its astrophysical significance in rotating black hole 
magnetospheres, and the basic formulae discussed in Sec. II for the stability analysis in Kerr geometry should be fully 
applied. The results of this paper would be very useful as the foundation of the future problems. 
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